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An averaging method is applied to derive effective approximation 
to the following singularly perturbed nonlinear stochastic damped 
wave equation 

uutt + ut = Au + f{u) + v'^W 



on an open bounded domain D C M" , 1 < n < 3 . Here z^ > is 
a small parameter characterising the singular perturbation, and v'^ , 
< Q < 1/2 , parametrises the strength of the noise. Some scaling 
^ ■ transformations and the martingale representation theorem yield the 

'^ . following effective approximation for small v, 



00 



f— ^ ■ to an error of o(z^") . 



ut = Au + f{u) + i^^W 
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1 Introduction 

Wave motion is one of the most commonly observed physical phenomena, 
and typically described by hyperbolic partial differential equations. Non- 
linear wave equations also have been studied a great deal in many modern 
problems such as sonic booms, bottlenecks in traffic flows, nonlinear optics 
and quantum fleld theory [131 [201 e.g.]. However, for many problems, such as 
wave propagation through the atmosphere or the ocean, the presence of tur- 
bulence causes random fluctuations. More realistic models must account for 
such random fluctuations. Hence we study stochastic wave equations (TJ [U 
e.g.]. 

Here we study an effective approximation, in the sense of distribution, for 
the following nonlinear wave stochastic partial differential equation (spde). 
The SPDE is a singularly perturbed problem on a bounded open domain 
D cW\l<n<3, 

z/<i + < = AM'^ + /(M'^) + z/°iy, m'^(0) = uo, <(0) = ui, (1) 

with zero Dirichlet boundary on D. Here z/° with < z/ < 1 and < a < 1/2 
parametrises the strength of noise, and A is the Laplace operator in M". The 
noise W{t) is an inflnite dimensional Q- Wiener process which is detailed in 
section [2j The spde ([T]) also describes the motion of a small particle with 
mass u and an infinite number of degrees of freedom [U [5] . We are concerned 
with the effective approximation of the solution to the spde (JT]) for small 
z/ >0. 

For a = 1/2 , the limit of the random dynamics of spde ([T]) as z/ — )■ has 
been studied by Lv and Wang [TTl p^8]. The random attractor and measure 
attractor of spde ([T]) are approximated by those of the deterministic pde 

Ut = Au + f{u) (2) 

as z/ — )■ in the almost sure sense [TT] and weak topology [IB] respectively. 

The important case of a = , which is an infinite dimensional version of 
the Smolukowski-Kramers approximation, is studied by analysing the struc- 
ture of solution of linear stochastic wave equations [U |5]. For any T > , 
the solution u{t) to the spde ([T]) is approximated in probability by that of 
the stochastic system 

ut = Au + f{u) + W 

as z/ ^ in space C(0, T; L'^{D)). 

Here we extend the approximating result to the case when < a < 1/2 
and derive a higher order approximation in the sense of distribution. Re- 
cently, the stochastic averaging approach was developed to study the effective 



approximation to slow-fast SPDEs [HI H] in the following form 

V \ V 



where / and g are nonlinear terms, 0\ and 02 are some constants, and 
W\ and W2 are Wiener processes. Notice that upon introducing v" = u^, the 
SPDE ([T]) is rewritten as 

v!; = v", m'^(O) = uq , 

which are also in the form of slow-fast SPDEs. Then we can follow the stochas- 
tic averaging approach to derive an effective averaging approximation of u'^ , 
the solution of SPDE ([T]) as i/ — )■ for all < a < 1/2 . Here the case 
a = 1/2 is the most important case because all cases of a G [0, 1/2] can 
be transformed to the case a = 1/2 , see section H] and section O By an 
averaging approach and martingale representation theorem, we prove that 
for small z/ > with < a < 1/2 the solution of SPDE ([T]) is approximated 
in the sense of distribution by u'^ which solves 

u'^ = Au" + fiu") + u'^W , m'^(O) = Mo , (3) 

where W{t) is a Wiener process distributes same as W{t). This result shows 
that for any small z/ > with < a < 1/2 the term i'u'({t) is a higher order 
term than the random force term v°'W{t). 

Section [3] gives the approximation for the important case that a = 1/2 . 
Previous research ^JJ gives an approximation which is a deterministic equa- 
tion. However, our new approximation shows that for small z/ 7^ , there is 
a small fluctuation which distributes same as y/uWit), see (j3]). This gives a 
more effective approximation. 

Section[6]explores a parameter regime where a nonlinear coordinate trans- 
formation underlies the existence of a stochastic slow manifold for the case 
a = 0. The stochastic slow manifolds of both the SPDE ([T]) and the model ([3]) 
have the same evolution in the parameter regime and so provide evidence of 
the stronger result of pathwise approximation therein. 



2 Preliminaries 

Let D C M" , 1 < n < 3 , be a regular domain with boundary F. Denote 
by LF'{D) the Lebesgue space of square integrable real valued functions on D, 



which is a Hilbert space with inner product 

{u,v) = / u{x)v{x) dx , u,v G L'^{D) . 



J D 

Write the norm on L'^{D) by HmHo = {u,u)^^'^ . Define the following abstract 
operator 

Au = -Au , u e Dom(A) = {u e L^{D) : Au e L\D) , u\r = 0} . 

Denote by {A^} the eigenvalues of A with < Ai < A2 < ■ • • < A^ < ■ ■ ■ , 
Afc — )■ 00 as A; — 7- 00 . For any s e M , introduce the space Hq{D) = Dom(A'^/^) 
endowed with the norm 

\\u\\s=\\A'/\\\o, ueH^{D). 

Consider the following singularly perturbed stochastic wave equation with 
cubic nonlinearity on D: 

z/< + < = AM'' + /3M"-(M")3 + z/"Vr(t), (4) 

u^iO) = uo, <(0) = Mi, (5) 

u'\r = 0, (6) 

with < 1/ < 1 and < a < 1/2 . Here {W{t)}tm is an L'^ {D)-valued two 
sided Wiener process defined on a complete probability space (fi, J^, {J-'t}t>o, P 
with covariance operator Q such that 

Qek = bkCk, k = 1,2,... , 

where {ck} is a complete orthonormal system in L'^{D), b^ is a bounded 
sequence of non-negative real numbers. Then 

00 
W{t) = ^^/hekWkit), 
fc=i 

where Wk are real mutually independent Brownian motions [T2]. Further, we 

assume 

00 00 

Bo = '^bk < 00 and i?i = ^ Xkbk < cx) . (7) 

k=l k=l 

Then by a standard method [8j, for any {uo,Ui) G H^'^^{D) x ^{D), s G M, 
there is a unique solution u^ to ©-(E]), 

u^eL\n,C{Q,T-H^,-^\D))), (8) 

<GL2(fi,C(0,T;i7^(D))). (9) 

In the following we write f{u) = (5u — v? and F{u) = J^ f{r) dr . 
For our purpose we need the following lemma. 



Lemma 1 (Simon [17]). Assume E, Eq and Ei be Banach spaces such that 
El (& Eq , the interpolation space {Eq, -Ei)6»,i C E with 6 G (0, 1) and E C 
Eq with C and ^ denoting continuous and compact embedding respectively. 
Suppose po,Pi G [1, oo] and T > , such that 

V is a bounded set in L^^(0, T; Ei), and 

dV := {dv.ve V} IS a bounded set m LP°{0,T; Eq). 

Here d denotes the distributional derivative. If 1 — 9 > 1/pe with 

1 _ 1-^ 6 
Pe Po Pi ' 
then V is relatively compact in C{0,T; E). 

In the following, for any T > , we denote by Ct a generic positive 
constant which is independent of u. 



3 The case of a = 1/2 



We first consider the special case of a = 1/2 which was recently studied by a 
direct approximation method fUX [T8] . Here we apply an averaging method 
to give more effective approximation to equation (JH)-©. We rewrite dl])-© 
in the form of slow-fast SPDEs: 



du" 
dv" 



-- [v'' - Au" - /K)] dt + ^ dW{t) 

U \ v 



V\^) = Ml . 



Notice that the slow part u'^ and fast part v'^ are linearly coupled, 
simplicity we consider {uo,Ui) e {H'^{D) n H^{D)) x H^{D). Then (g 
has a unique solution in L\n,C{0,T; (H^D) f] H^{D)) x H\D))). 



(10) 

(11) 

For 



3.1 Tightness of solutions 

Let (m^, v'^) be a solution to flTOl) - flTTl) with z/ > . In order to pass to the 
limit z/ — )■ in the averaging approach, we need some a priori estimates on 
the solutions. 

Theorem 1. Assume Bi < oo . For any T > 0, there is a positive con- 
stant Ct such that 



E 



max ||M''(t)||2 + max ||f''(t)||o 



0<t<T 



0<t<T 



<Ci 



:i2) 



and for any integer m > 



E [ \\u''{t)\\l"'dt<CT. 
Jo 



Proof. The result on ||'u'^(t)||2 is found by a simple energy estimate [TS]. Now 
we give the estimate on ||t>'^(t)||o. By equation (ITT]) . 

v''{t) = e-'/^ui + - / e-^'"^/" [Am"(s) + /K(s))] ds 
^ Jo 

+ ^ / e-^'-'^/^dWis). 

'y .In 



Noticing assumption ([7]), by the estimate on ||M'^(t)||2 and maximal inequality 
of stochastic convolution p^ Lemma 7.2], 



E 



0<t<T 



max II t> i(.;i|o 



<C7 



for some positive constant Ct- The last inequality of the theorem is obtained 
by the same method [IB] and Poincare inequality. This completes the proof. 

D 

Now by the above estimates and LemmalU we have the following theorem. 

Theorem 2. For any T > 0, {C{u'^)}q^^<i the distribution of u^ is tight in 
the space C{0,T;H^{D)). 

By the above tightness result, to determine the limit of u'^ we can pass 
to the limit i/ — )• in a weak sense; that is, for any ip G C^{D), we consider 
the limit of M^''^(t) = (u^it), if) in the space C(0, T) as z/ -> . 

3.2 Limit of m^'^ in C(0, T) 

Now we pass to the limit i/ — ;■ in {u'^''^} in the space C(0, T) for any T > . 
First, by equations (ITU1) - (ITT1) . 

du"'^ = v''''^dt, (13) 

dv"''' = --K■^ + (Vn^V^)-(/K),<^)]rft+^rfly'^(t), (14) 
u \ v 



with M^''^(0) = (mo, <^) and t;^''^(0) = (mi, y?) where v""^^ = {v", ^) and W^{t) = 
{W{t),ip). In the following we also write v'^''^ as t;'^''^'"W which shows the 
dependence of v'^''^ on the slow part u'^. 

6 



Second, for any fixed u G H'^{D) fl Hq{D) we consider the fast equation 
dv"'"" = -- [v"''' - An - /(n)] dt + ^ dW{t) . (15) 



Equation flT5|) lias a unique stationary solution v^'^ . Moreover, the stationary 
solution v^'"^ is exponentially mixing and the distribution of v"'"^ is the normal 
distribution M {^u^ f{u),QI1) [4]. 

Now for any u € }i'^{D) n lil{D) define 

/7^(u, t) = z/ [t;'^'"(t) - t;^'"(0)] + /" K'"(s) - An - /(n)] ds . 
Then n"^'*^ solves the following equation 

n^'^(t) = (no, ^) - y [(Vn^(s), Vv^) - (/(n^(s)), ^)] ds 

+ (/J'^(n'^(t), t), y.) - z/ (t;'^'"(t) - ^'''"(O), y.) . (16) 

Third, we study the behaviour of {H'^{u'^{t),t), ip) for small v. Let H'^''^{u, t) = 
{H^^u, t),ip) , then define 

M^''' = ^H''''^{u''{t),t). (17) 



By the definition of H'^''^{u,t) and equation flTSl) . M^'^''^ is a martingale with 
respect to {J-^ : t > 0}, and the quadratic covariance is {M^''^)t = t{Q(p,(p). 
Now define R'''^{t) = - (n^'"(t) - n'''"(0), (^), then rewrite (^ as 

n^'^(t) = (no,^)- f\{Vu''{s),Vy^) - (/(n^(s)), ^)] rfs + ^JM^ + z//?'^''^(t) . 

(18) 
Invoking Theorem [H 

0. (19) 



limE 



max v^|i?^''^(t)| 

0<t<T 



Then define the process 



Mr = ^\u'''''{t)-{uo,y^)+ I [{Vu%s),Vy^)-{f{u''{s)),^)]ds\. 



(20) 



By the definition of H'^''^{u,t) and fITI?]) we liave tlie tiglitness of A^^'"^ in 
space C(0, T) for any T > . Let P be a limit point of the family of prob- 
ability measures {C{J^^''^)}o^u<i and denote by Aif, a C(0, T)-valued ran- 
dom variable with distribution P. Let \l/ be a continuous bounded function 
on C{0,T). Set "^"{s) = ^(m'^'^(s)), then noticing (^, 

which yields that the process {A^f }o<t<r is a P-martingale with respect to 
the Borel a-filter of C(0,T). 

We consider the quadratic covariation of the martingale Aif. By the 
definition of A^^'*^, passing to the limit z/ — )• in f l2U]) . we derive A^f is a 
square integrable martingale with the associated quadratic covariation pro- 
cess is {Q(p,(p)t. Then by the representation theorem for martingales [TU] . 
without changing the distributions of Al^'"^ and Aif, one extends the origi- 
nal probability space {Q, J^, P) and chooses a new Wiener process Vr'^(t) such 
that Alf = a/QW^'^('^)) which is unique in the sense of distribution. 

By the definition of Al^''^, VT*^ can be chosen as {W, (p) where VT is a 
cylindrical Wiener process. Then from fl20p we have in the sense of distribu- 
tion 

K(t),^) 
= (wo, ^)- I [(Vw^(s), V^) - (/K(s)), ^)]ds + v^Air + o(v/^) 

for any ip G C^{D). Then by discarding the higher order term and the 
tightness of u^ , we have the following approximating equation 

dv!" =[Au'' + f{v!')]dt + ^dW^ , (21) 

where W^ is some an L'^{D) valued Q- Wiener process. 

Theorem 3. Assume Bi < oo and a = 1/2 . For small u > , there 
is a new probability space {Q,^,F), an extension of the original probabil- 
ity space {Q,J^,F), such that for any T > 0, the solution u" to / fJOj) - 
ni\) is approximated by u^ which solves l[2l\). to an error of o[^\jV) , in the 
space C(0, T; Hq{D)) for almost all u E fl . 

The above spde (pTj) is more effective than the limit pde (j2]) [TT] as it 
incorporates fluctuations for small z/ > . This result also implies that the 
singular term I'u^it) is a higher order term than y/vW{t) for small z/ > in 
the sense of distribution at least. The following two sections show that z/ti^(t) 
is always a higher order term than u'^Wit) for any < a < 1/2 . 



4 The case of a = 

Next we consider the case of a = ; that is, consider the following spde 



z/<, + < = Au'^ + Pu'' -{u^f + Wit), 

u'-'iO) = Mo, <(0)=Mi, 

0. 



«1r 



(22) 
(23) 
(24) 



First we have the following a priori estimates on u'^ in the space C(0, T; Hq{D)). 

Theorem 4 (Cerrai & Freidlin [5])- Assume Bi < oo . For any T > , there 
is a positive constant Ct such that 



E 



max \\u 

0<t<T 



<Ci 



We follow the approach for the case of a = 1/2 . For this we introduce 
the scalings u'^ = \fvu^ and v^ = \fi'u\ . Then 



dv" 



V^dt, M^(O) = 0^Mo, 



1 



r - Am" - v^/ 



u 



dt + ^dW{t), v''{0) = ^/^Ul. 



V 



By standard energy estimates [18], by a similar discussion to that in Section [3], 
and by Theorem HI we have the following theorem. 

Theorem 5. Assume Bi < oo . For any T > , there is a positive con- 
stant Ct such that 



E 



max ||M'^(t)||9+ max llf^ 
o<t<r o<t<T 



<Ci 



and for any integer m > 



E 



[ \\u''(t)\\l"'dt<CT. 

Jo 



Moreover, the distribution of u'^ is tight in space C{0,T] Hq{D)). 

We consider the asymptotic approximation of u" for small z/ > . For 

any cp G 0^(0), let u"'^ = (u^^if), v"''^ = {d\ip) and W^it) = {W{t),ip). 
Then 



du"'^ 
dv"''^ 



d^'^dt 
1 



[m"'^ + (VM^ V^) - v^(/(m7v/^), ^)] dt+— dW^{t), 



with ^'^•^(0) = (m^(0),(^) and 'D^'^(O) = (v^iO),^). 

We also consider the following fast spde for fixed u and u G H'^{D) fl 

dv"'^ = -- [v"'^ - Am - v^/ (m/v^)] dt + ^ dW^(t) . (25) 

U \ v 



For fixed z/ G (0, 1] and u G H'^{D)r\H^{D), SPDE ([22]) has a unique station- 
ary solution with the normal distribution M {/S.u + y/D f {u / y/v) , Q/2) |1]. 
Now for any u G H'^{D) n i/o(£') define 

H^iu.t) = V [{}'^'"(t) - {}"'"(0)] + / [{}"'"(s) - Am - v^/(m/v^)] t/s . 

Thus we can follow the same discussion in last section for the case of a = 1/2 . 
We write 

S-'^(t) = y^(„o, ^)- I {Vu\s), V^)ds + V^ [ {f {u%s)/V^) , v)ds 

Jo Jo 

+ V^Mr, (26) 

where ^/uAit''^ is the remainder term. By a similar discussion to that of the 
last section, Jli'^''^ is tight in space C(0,T) for any T > . Let P be a limit 
point of the family of probability measures £{A^^''^}o<!/<i in space C(0,T). 
Let Aif be a C(0,T)-valued random variable with distribution P. Then we 
have the following lemma. 

Lemma 2. For any ip G C^{D), the process Aif defined on the probabil- 
ity space {C{0,T),B{C{0,T)),P) is a square integrable martingale with the 
associated quadratic covariation process {Qf, (f)t . 

By the representation theorem for martingales (TU], without changing 
the distributions of A^)|'''^ and Aif one can extend the original probability 
space {Q,J^,F) and choose a new cylindrical Wiener process W{t) such that 
J\4f = ■\/Q{W, (f), which is unique in the sense of distribution. 

Then in the sense of distribution by (l26l) we write out 

(M^(t),y;) = ^{u^,ip)- f {Vu''{s),Vip)ds + V^ f {f{u''{s)l^),ip)ds 

Jo Jo 

= v^(mo, <^) - / (Vn'^(s), V^)ds + V^ [ {f (u"(s)/v^) , y^)ds 
Jo Jo 

+ V^VQ{W,^)+o{V^) (27) 

10 



for any yj G C^{D). Then we have, noticing that m*^ = y/uu", the following 
approximating SPDE for small u > 0: 

du'' = [Au" + f{u'')]dt + dW^, u^iO) = Uq , (28) 

where W^ is some L'^{D) valued Q- Wiener process. Then we infer the fol- 
lowing result. 

Theorem 6. Assume Bi < oo and a = . Then for small u > , there 
is a new probability space {Q, ^, P) which is an extension of the original 
probability space {Q, J-", P) such that for any T > , the solution u" to / f^) - 
([^P is approximated by u" which solves ^EB\), to an error of o{l), in the 
space C(0, T; Hq{D)) for almost all u E Cl . 

5 The case of < a < 1/2 

Now we consider the case of < a < 1/2 ; that is, consider the following 

SPDE 

z/<, + < = Au'' + (3u''-{u''f + u''W{t), (29) 

u^iO) = no, <(0) = ni, (30) 

u^lr = 0. (31) 

First, by the same analysis as Theorem IU we also have the following 
result on the a priori estimates on u'^. 

Theorem 7 (Cerrai & Freidlin [S]). Assume Bi < oo . For any T > , there 
is a positive constant Ct such that 



E 



max IIm'''''^^"^ 

0<«T 



<Ct. 



We also apply the method in Section [31 Make the following scaling trans- 
formation n!^ = u^/'^~°'u'^ and v''' = v^l'^~'^v^ . Then 



du^ = v^dt . 
dv" = -- 

V 






dt+^dW(t), 

V 



By a direct energy estimate or the scaling transformation and Theorem [7] we 
deduce the following theorem. 

11 



Theorem 8. Assume Bi < oo . For any T > , there is a positive con- 
stant Ct such that 



E 
and for any integer m > 



max \\u''(t)\\l+ max ||'y''(t)||n 

0<t<T ^ 0<t<T ^ 



<Ct, 



E 



/ \\u''{t)\\l'^dt<CT. 

Jo 



Moreover, the distribution of u'^ is tight in space C{0,T; Hq{D)). 

Then we can follow the same discussion of Section [Hand have the following 
result. 

Theorem 9. Assume Bi < oo and < a < 1/2 . For small u > , there 
is a new probability space {fi, ^, P) which is an extension of the original 
probability space {Q, J-", P) such that for any T > , the solution u'^ to / f^) - 
(EIP is approximated by u'^ which solves 



du" = [Au^ + f{u^)]dt + u'^dW^, u^{0) = uo , (32) 

to an error of o{h'°'^ , in the space C{0,T; Hq{D)) for almost all tu E 0.. 

6 A stochastic slow manifold compares the 
SPDEs for the case of a = 

This section shows the long time effectiveness of the averaged model by com- 
paring it to the original via their stochastic slow manifolds. 

We compare the spde (1221) and its model spde fl28p in a parameter 
regime where both have an accessible stochastic slow manifold. Consider 
the SPDE fl22l) restricted to one spatial dimension as 



uutt + ut = u^x + fiu) + aW where / = (1 + l3')u - u^. (33) 

Consider this spde on the non-dimensional domain D = (0, vr) with bound- 
ary conditions M = 0onx = 0,7r. The parameter a here explicitly measures 
the overall size of the Q- Wiener process W{t) which by (jTj) is finite. The 
small parameter /3' measures the distance from the stochastic bifurcation 
that occurs near /3' = . In this domain there will be a stochastic slow man- 
ifold of the SPDE fl33l) that matches the slow dynamics in the approximating 
SPDE (j28j) . This section compares the stochastic slow manifolds. 
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The SPDE fl5^ has a technically challenging spectrum. However, the con- 
struction of its stochastic slow manifold is easiest by embedding the spde fl55]) 
as the 7 = 1 case of the following slow-fast system of spdes 

Ut = U:,::,: + U + V , (34) 

uvt= -v--fiy {d^^ + 1) Mj + /3'u -u^ + aW . (35) 

The parameter 7 controls the homotopy: from a tractable base when 7 = 
as then all linear modes in the very fast v equation (135|) decay at the same 
rate l/u (and the slow u modes of sin /ex have decay rates 1 — k"^); to the 
original SPDE 033 p when 7 = 1 (upon eliminating v). 

A stochastic slow manifold appears On the non-dimensional inter- 
val (0,7r), with Dirichlet boundary conditions on u, the eigenmodes must 
be proportional to sin kx for integer wavenumber k. Neglecting noise tem- 
porarily, cr = in this sentence, for all u < 1 and all homotopy parame- 
ter < 7 < 1 there is one zero eigenvalue and all the rest of the eigenvalues 
have negative real part; the slow subspace corresponding to the neutral mode 
is spanned by {u,v) oc (sinx,0) (local in {u,v,a), but global in u and 7). 
By stochastic center manifold theory [H [3] , and supported by stochastic nor- 
mal form transformations [2l[T5l|16], when the noise spectrum truncates and 
the nonlinearity is small enough, the dynamics of the spdes (!34|) - (!35|) are 
essentially finite dimensional and a stochastic slow manifold exists which is 
exponentially quickly attractive to all nearby trajectories. 

Computer algebra constructs the stochastic slow manifold We seek 
the stochastic slow manifold as a systematic perturbation of the slow sub- 
space u = a sin x . The intricate algebra necessary to handle the multitude of 
nonlinear noise interactions is best left to a computer [Hi [161 ^.g.]. However, 
the following expressions may be checked by substituting into the govern- 
ing SPDES (jMD^dSSD and confirming the order of the residuals is as small as 
quoted — albeit tedious, this check is considerably easier than the derivation. 
The evolution on the stochastic slow manifold may be written 

a = (3'a- |a^ + [l - 2u(3' + \va^ - -^^a^\ hwi 

(36) 

The stochastic slow manifold itself involves Ornstein-Uhlenbeck processes 
written as convolutions over the past history of the noise processes: define 
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e f^^-kw = J_^exp[—fi(t — s)]dws for decay rates /i^ = /c^ — 1 characteristic 
of the fcth mode. Then the stochastic slow manifold is 

n = a sin X + -^a^ sin 3x — -^a^ [b^e^^'^-kws sin x + bie^^^-kwi sin 3a;] 
+ y^ bk [l + fik^ + -yi^ifik - Aifce"'^'=V)] e~^''^-kWk sin kx 

k>2 

- ^ bkC^^^^i^Wk sin fcx + /3' ^ bke'^^^i^e^^^^i^Wk sin fcx 

A;>1 fc>2 

+ 4 5Z {^fc+2e"'"'**e"^'=+2Vwfc+2 sin kx - 2bke~^''^i<e~^''^-kWk sin kx 

k>2 

+ 6fce-'^'=+2Ve-'^'=**Wfc sin[(A; + 2)x] } + 0{u^ + f3'^ + a\ a^) , (37) 

and a correspondingly complicated expression for the field v{x,t). Observe 
that the slow SDE (I5B]) does not contain any fast time convolutions from 
the Ornstein-Uhlenbeck processes: it would be incongruous to have such 
fast processes in a supposedly slow model. We keep fast time convolutions 
out of the slow SDE f l36|l by introducing carefully crafted terms in the slow 
mode sinx in the parametrization of the stochastic slow manifold fl37|) : here 



the amplitude of the slow mode sin a; is approximately a — ^a'^b^e'^^-kw^ — 
bie^^^'^'kwi . Other methods which do not adjust the slow mode either average 
over such adjustments and so are weak models, or invoke fast processes in 
the slow model. 

Note that the homotopy parameter 7 affects the stochastic slow manifold 
shape (1371) . but only weakly. To this order the homotopy has no effect on 
the evolution on the stochastic slow manifold 



Compare with SPDE (125]) The corresponding stochastic slow manifold 
of the SPDE fl28p . in this parameter regime, is straightforward to construct, 
via the web server pjij for example. For stochastic slow manifold u ~ a sin a; 
one finds the corresponding slow SDE 

a = (3'a- |a^ + [l - i^*^^] ^i""^! 

+ M+ lis/^')^' - lii^l ^3^3 + Tk^%w, + o(/3" + a^ a) . (38) 



This slow SDE is symbolically identical with the SDE fl36l) . one just removes 
the overbars. We conclude that these stochastic slow manifolds confirm the 
modeling of the SPDE fl^ by its model SPDE (1251) : at least in the regime of 
one space dimension with small amplitude a, bifurcation parameter /?', and 
finite truncation to the noise. 



14 



Acknowledgements The research was supported by the NSF of China 
grant No. 10901083, Zijin star of Nanjing University of Science and Technol- 
ogy, and the ARC grant DP0988738. 

References 

[1] Ludwig Arnold. Random Dynamical Systems. Springer Monographs in 
Mathematics. Springer, June 2003. 

[2] Ludwig Arnold and Peter Inikeller. Normal forms for stochastic dif- 
ferential equations. Probab. Theory Relat. Fields, 110:559-588, 1998. 
ht tp : //dx ■ doi . org/10 . 1007/s00440005015 9 , 

[3] P. Boxler. A stochastic version of the centre manifold theorem. Probab. 
Th. Rel. Fields, 83:509-545, 1989. 

[4] S. Cerrai & M. Freidlin, On the Smoluchowski-Kramers approximation 
for a system with an infinite number of degrees of freedom, Prob. Th. 
and Relat. Fields 135 (2006), 363-394. 

[5] S. Cerrai & M. Freidlin, Smoluchowski-Kramers approximation for a 
general class of SPDEs, J. Evol. Equa. 6 (2006), 657-689. 

[6] S. Cerrai and M. Freidlin, Averaging principle for a class of stochastic 
reaction-diffusion equations, Probab. Theory Relat. Fields 144 (2009), 
137-177. 

[7] P. L. Chow, Stochastic wave equation with polynomial nonlinearity, 
Ann. ofAppl. Prob. 12(1) (2002), 361-381. 

[8] P. L. Chow, Stochastic Partial Differential Equations. Chapman & 
Hall/CRC, New York, 2007. 

[9] P. L. Chow, W. Kohler, & G. Papanicolaou, Multiple Scattering and 
Waves in Random Media, North-Holland, Amsterdam, 1981. 

[10] N. Ikeda and S. Watanabe, Stochastic Differential Equations and Diffu- 
sion Processes, Vol. 24 of North-Holland Mathematical Library, North- 
Holland Publishing Co., Amsterdam, 1981. 

[11] Y. Lv & W. Wang, Limit dynamics for stochastic wave equations, J. 
Diff. Equa. 244 (2008), 1-23. 



15 



[12] G. Da Prato & J. Zabczyk, Stochastic Equations in Infinite Dimensions, 
Cambridge University Press, 1992. 

[13] M. Reed & B. Simon, Methods of Modern Mathematical Physics II, Aca- 
demic Press, New York, 1975. 

[14] A. J. Roberts. Resolving the multitude of microscale in- 

teractions accurately models stochastic partial differential equa- 
tions. IMS J. Computation and Maths, 9:193-221, 2006. 



http : //www ■ 1ms . ac . uk/ j ciii/9/lms2005-032 



[15] A. J. Roberts. Normal form transforms separate slow and fast modes in 
stochastic dynamical systems. Physica A, 387:12-38, 2008. 

[16] A. J. Roberts. Normal form of stochastic or determin- 

istic multiscale differential equations. Technical report. 



http : //www . maths . adelaide . edu . au/anthony . roberts/sdenf . html 



2009. Revised April 2011. 

[17] J. Simon, Compact sets in the space L^(0,T; B), Ann. Mat. Pura AppL, 
146 (1987), 65-96. 

[18] W. Wang and Y. Lv, Limit behavior of nonlinear stochastic wave equa- 
tions with singular perturbation. Disc, and Cont. Dyna. Syst. B, 13(1) 
(2010) 175-193. 

[19] W. Wang and A. J. Roberts, Average and deviation for slow-fast 
stochastic partial differential equations, preprinted, 2009. 

[20] G. Whitham, Linear and Nonlinear Waves, Wiley, New York, 1974. 



16 



